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<D ■ Abstract 



In this paper, we address the possibility of generalising the standard analysis of thermal contact 
between a sample system and a heat bath, by including long range interactions between them. 
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As a concrete example, both system and bath are treated within the long range Ising model. For 

a ' 

■ this model, we derive the equilibrium probability distribution of the energy of the sample system. 

Equilibrium properties of the system magnetisation and stability of the solutions is discussed. We 

o , 

O , find existence of a metastable phase below a critical temperature of the bath. 
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I. INTRODUCTION 



The theory of equilibrium thermodynamics and statistical mechanics has been rigorously 
developed only for the case of short range interactions among the components of the system. 
The presence of long range interactions among different components, makes a system much 
more complex and the standard approaches become inapplicable. Currently, there is a lot 
of interest to develop new methods and tools to deal with systems involving long range 
interactions 111. Apart from the standard examples within the fields of cosmology and 
astrophysics [2j, a growing number of physical laboratory systems have recently emerged in 
which the interactions are long-ranged, notably in the areas of plasma physics nuclear 
physics and atomic clusters jjj, Bose-Einstein condensates jjj, 2D hydrodynamics 0, Q], 
magnetism and so on (see [l| and references therein). In systems with long 

range interactions, energy is generally nonadditive. Due to this feature, the thermodynamics 
of such systems displays unusual properties, like inequivalence of different ensembles 12. 1^|. 



possible temperature discontinuity at first order transitions and negative specific heat in 
microcanonical ensemble [l^j ]. 

Now the standard derivation of Boltzmann-Gibbs distribution considers a thermal contact 
between a sample system and a heat bath with a very large heat capacity. The derivation 
assumes an additive property for energy, to arrive at the desired result. The case of long range 
interactions between the system and bath seems to be beyond the scope of this derivation 
Q]. For one, to include such interactions, we have to specify the nature of the interactions 
between the components of the system and the bath. This implies specifying the microscopic 
structure of the bath also which makes the characterisation of the bath more detailed. 
In contrast, within canonical ensemble the bath is characterised only by its temperature 
irrespective of its microscopic model. Further, nonadditivity of energy brings additional 
difficulties. 

In this paper, we consider this latter case and derive the probability distribution for a 
system in thermal contact with a heat bath in the presence of long range interactions. In 
section II, we take a specific model for long range interactions of mean-field type between 
the system and the bath, and present such a derivation. In section III, we compare our 
model with the Curie- Weiss model which is paradigmatic model for ferromagnetism. In 
section IV, equilibrium value of the magnetisation is discussed and stability of the solutions 
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is highlighted. The last section V, presents a summary of our approach and results. 



II A. Canonical Ensemble from Microcanonical Ensemble 

For the sake of completeness and to fix the notation, let us review one of the text-book 
derivations of the canonical ensemble [l£]. Denote by E\ the energy of the sample and by E 2 
the energy of the reservoir. The sample and the reservoir together form an isolated system 
and the interactions between them are usually considered to be short ranged. Thereby the 
total energy of the composite system say E , is given by the sum of energies E% + E 2 and 
is constant. Next, fli(Ei) where i — 1,2, denote the respective number of microstates at 
the given values of energies. The probability that the system 1 is in a certain microstate of 
energy E\ is given by 

where E 2 = Eq — E±, and Qi +2 (Eq) is the total number of states available for the composite 
system 1 + 2. Thermodynamic entropy of the reservoir is given by Boltzmann's formula: 

S 2 (E 2 ) = k B \nn 2 (E 2 ). (2) 

Therefore, pi(E 1 ) oc fii(-Ei) e 52//feB . In the following, we keep by k B = 1 for simplicity. Now 
expand S 2 (E 2 ) around the most probable value E 2 as: 

S 2 (E 2 ) = S 2 (E 2 ) + (^j _ (E 2 - E 2 ) + ■ ■ ■ . (3) 

As the assumed size of the bath is quite large, its instantaneous energy will be very close to 
the most probable energy. Therefore, taking E 2 —E 2 to be small, the higher order terms in the 
above expansion are negelected. Let the inverse temperature parameter {dS 2 / dE 2 ) ^ 2 = (3 2 . 
Then using additive property of the energy, we have E 2 — E 2 = E\ — E\ and thus we obtain 
the Boltzmann distribution for system energies: p{Ei) oc Qi(Ei) exp(—/3 2 Ei). 

IIB. ISING MODEL WITH LONG RANGE INTERACTIONS 

Now consider the Ising spin model for both the system and the bath interacting with 
each other through long range interactions. (For a recent discussion of mutual thermal 



equilibrium and alternate thermodynamic descriptions within this model, see Ref. Q-) 
Take a lattice of iV spins with the hamiltonian given by 

E = -Y. h 3 a 3- E fjk a J ( 4 ) 

3 <jk> 

Each spin variable cr, is defined over the set {+1,-1}. The total system is assumed to be 
separated into two regions representing a sample system and a heat bath. The spin excess 
in each region is given by = — N[~ , where are the number of up (+) and down 
(-) respectively. Jjk and hj are the known coupling constants within each region which for 
simplicity, are taken to be constant over a region. We focus on the ferromagnetic interactions 
which mean J7"'s are greater than zero. In terms of new variables Ei = —hiCi, sometimes 
called the local energies, the total energy upto a constant term, can be written in the form 

E = E x + E 2 - aEx 2 - bE 2 2 - cE x E 2 , (5) 

where a, b and c are suitably defined constants to be elaborated later. Clearly, in the presence 
of long ranged interactions, the energies as well as the temperatures of the systems, get 
modified. For example, the most probable energy of system 1 is t\ = Ex — aEx and its 
inverse temperature is 

g* = ^l = §1 (a) 

Pl ~d£x (l-2aExY [) 

where fix = dSx / dEx is the inverse temperature in the absence of long range interactions 

within system 1. Similarly, for the bath we have 

P * 2 = (1 -2bE 2 ) ' (7) 
In the following, we adapt the scheme as given in Eqs. ©-(JB} to the case when the total 
energy is given by (|5|). It is convenient to express the entropy of a system in terms of the 
respective local energy E^. To evaluate (E 2 — E 2 ) for this case, we proceed as follows: 
In terms of the most probable values of the local energies, we have 



E = Ex + E 2 - aEx 2 ~ bE 2 - cE x E 2 . (8) 
Equating Eqs. © and ((HJ) and defining AEi = Ei — Ei, i = 1, 2, we can write 
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A£i + AE 2 - a(£i 2 - Ei) - b{E 2 + E 2 ).AE 2 - c{E x E 2 - E X E 2 ) = 0. (9) 

Now due to the large size of the system 2, the energy E 2 is very closely approximated by its 
most probable value. Thus the following terms in the above equation can be approximated 

as 

(E 2 + E 2 ).AE 2 ~ 2E 2 .AE 2 , (10) 

and 

{E X E 2 — EiE 2 ) ~ E 2 .AE 1 . (11) 
Using the above relations in Eq. ([9*|1. we get the expression 

= -A El (l -cE^ + ajES-ES) 

(l-2bE 2 ) ' 1 ; 

Finally, substituting f)12|) in (jSJ), we arrive at the probability distribution for energies of 
system 1, 

P(E 1 ) = ' h[hl) , (13) 

where 

Z = £ Sl x {E x ) exp (-${(1 - cE 2 )E 1 - aE 1 2 }) , (14) 

Ei 

the sum being over the range of energies accessible to system 1. It is appropriate here to 
explain the various paramters above: 

a = oi o— ™, c— — - — — . (15) 

2N 1 h 1 2 2N 2 h 2 2 N 2 h x h 2 y ' 

The scaling of parameters by the number of particles is in the same spirit as Kac's priscription 



17| and it guarantees that the effective hamiltonian (see Eq. (JTfij) below) is extensive and 
thus a proper thermodynamic limit is ensured. 
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III. ANALOGY WITH CURIE- WEISS MODEL 



Clearly, when there are no interactions with or within the bath, i.e. J 2 = J12 — implying 
b = c = 0, then we obtain the probability distribution of the Curie- Weiss model which 
is described within canonical ensemble with hamiltonain H = E\ — aE 2 and the bath at 
inverse temperature /3 2 - 

As mentioned above, the presence of long range interactions within the bath has the 
effect of modifying the bath temperature to /3 2 . On the other hand, the long range interac- 
tions between the system and the bath are incorporated effectively as a modified external 
magnetic field. This can be clearly seen by noting that the distribution (|13|) corresponds to 
an effective hamiltonian H' = (1 — cE 2 )Ei — aE\ . In terms of the magnetisation per spin, 
Hi = J^i&i/Ni = —Ei/{hiNi) and using definitions (|T5j) for the parameters, we can write 



H'( yi ) = -N x 



(hi + Jviv-i)y\ + \ J iVi 2 



(16) 



indicating the extensive property of the effective hamiltonian. Note that the effective applied 
field is (h\ + J12V2) = (1 — C-E^/ii, where h% is the actual applied field. We note that in 
contrast to the canonical ensemble, in the present case, the bath has a specific microscopic 
realisation. The expression for the entropy of bath given by Eq. (J2J), corresponding to the 
most probable value of magnetisation per spin y 2 is 



5*2(2/2) = ~N 2 



1 + y 2 ) ^ (1 + y 2 ) | (1 - y 2 ) ^ (1 - fe) 



(17) 



— — 2 

The most probable energy of the bath given by £2 = E 2 — bE 2 , can also be written as 
£2(2/2) = —N 2 (h 2 y2 + \ J%y%)- Then the inverse temperature of bath is 

dS 2 ( dS2\ 1 



1% 



dy 2 \dy 2 



' Inf^V (18) 



2(h 2 + J 2 y2) \ 1 ~y2, 

Note the limiting values for the bath temperature: 

(i) For y 2 — > 1, we have f3 2 — >• 00, which corresponds to bath temperature of zero degree 
Kelvin. 



6 



(ii) For yi — > 0, /3% — > 0, or in other words, the bath is most disordered at arbitrarily high 
temperatures. 

IV. EQUILIBRIUM PROPERTIES AND STABILITY 

It is straightforward to calculate the free energy A(yi) from the partition function by the 
standard methods of Hubbard- Stratonovich transformation or the saddle point approxima- 
tion. Thus we obtain 



The stationarity condition dA/dyi = yields a self-consistent equation for the equilibrium 
magnetisation per spin of system 1: 



Usually, the analogue of this equation in the Curie- Weiss model is analysed for hi = case 
to infer the existence of a critical temperature, above which the system is paramagnetic and 
below which it is ferromagnetic. In the ferromagnetic phase, two values of magnetisation 
±|£/i| are equally allowed, which actually reflects the symmetry of the hamiltonian. 

In the present case, even in the absence of external field (hi = 0), there is an effective 
magnetic field J122/2 due to long range interaction between system and the bath. Thus at high 
temperatures, we have a unique minimum of free energy at a non-zero value of magnetisation. 
Thus the sample is magnetised even at high temperatures. As the temperature is lowered, a 
metastable solution for magnetisation also appears alongwith the global equilibrium solution. 
The limit of metastability (defined by the inverse temperature @2 a ^ which the metastable 
state appears while lowering the temperature of the bath) can be calculated by finding the 
point of inflexion, where both first and second order derviatives of free energy vanish. The 
latter condition yields 



In [2 cosh{ / 9*(/i 1 + J 12 y 2 + J m )}] . 



(19) 




(20) 




(21) 



Using the above condition alongwith Eq. (|20jh we can show that 



1 



(22) 



yi = ±, 1 
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From this equation, it is clear that J\ > 1. The condition, f}\ J\ = 1 implies that 
yi = 0. But for hi = 0, this stationary solution is satisfied only for J 12 = 0, (see Eq. (J2U|) 
above). On the other hand, for J 12 ^ 0, the stationary solution is non-zero. Thus we see that 
the condition can be satisfied by y± ^ for an inverse temperture $2 which implies 
that the bath temperature is lower than the critical temperature for the usual para-Ferro 
transition which obeys (3^"^ J\ = 1. The actual temperature for limit of metastability may 
be calculated from Eqs. (I2U|) and (1211. 



V. SUMMARY 



We have introduced a new kind of ensemble to generalise the usual treatment of a ther- 
mal contact between system and heat bath by including long range interactions between 
them. The derivation is motivated by the standard derivation of the canonical ensemble 
from the microcanonical ensemble. However, the crucial difference is the lack of additivity 
of the energy. This does not yield exponential Boltzmann distributions. The form of the 
distributions is strongly dependent on the form of hamiltonian of the total sample plus bath 
system. Thus incorporating long range interactions makes the problem more involved as we 
have to specify a microscopic model for the bath as well as the nature of long range inter- 
actions between the system and the bath. In this paper, we have treated a very simple case 
when the interactions can be modelled by long range Ising model. An appropriate scaling of 
the interaction parameters with system or bath size helps to obtain a thermodynamic limit 
for the system properties. The observations for this analysis are that the usual para-ferro 
transition appears to be suppressed, with the presence of a net magnetisation for the system 
at high temperatures, even in the absence of an applied field. This happens because the long 
range interaction with the bath provides an effective magnetic field. Moreover, as the tem- 
perature is lowerd, there appears a metastable state. This happens at a temperature which 
is generally lower than the critical temperature of para-ferro transition in the Curie- Weiss 
model. It is hoped that the approach presented in this paper will motivate further studies 
with other model long-range interactions, such as slowly decaying interactions. This may fa- 
cilitate the characterisation of thermodynamic behaviour in systems when interactions with 
the bath cannot be neglected. 
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FIG. 1: Fig. 1: Free energy (Eq. (|19j) ) plotted against system magnetisation for different bath tem- 
peratures, parameterised by the values of y~2] increasing values imply decreasing bath temperature. 
The other paramters are set at hi = 0, hi = 0.5, J\ = 4.0, J2 = 3.5, J\i = 0.05. 
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